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Abstract 

The A^{L) statistic characterizes the fluctuations of the number of levels as a function of the 
length of the spectral interval. It is studied as a possible tool to indicate the regular or chaotic 
nature of the underlying dynamics, to detect missing levels and the mixing of sequences of levels 
of different symmetry, particularly in neutron resonance data. The relation between the ensemble 
average and the average over different fragments of a given realization of spectra is considered. 
A useful expression for the variance of A^{L) which accounts for finite sample size is discussed. 
An analysis of neutron resonance data presents the results consistent with a maximum likelihood 
method applied to the level spacing distribution. 
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I. INTRODUCTION 



Neutron resonances provided the first context for modeling 



Dliysical reality with Random 



Matrix Theory (RMT) for a brief history of RMT see Bohr's compound nucleus 
description [sl identified the positions of the resonances as eigenvalues of the unknown com- 
plicated Hamiltonian governing the compound nucleus. Calculating the energies of these 
excited states is impossible, even for non-interacting particles the level density is prohibitive 
just from combinatorics alone; with interactions included, exact calculations even in a rea- 
sonably truncated space quickly become impossible. However, robust statistical features of 
the spectra are calculable. Statistical spectroscopy had already been opened up by Gurevich 
in 1939 when he investigated the regularities of level spacings in nuclear spectra ^. Wigner 
took it a step farther in 1951 when he suggested that although the specific energies cannot be 
calculated, the wave functions are so complicated that the statistical behavior of the energy 
levels mimics that of the eigenvalues of an ensemble of random matrices whose elements 
have probability distributions that do not favor any particular basis, i.e. their statistics are 
invariant under orthogonal transformations (change of basis). The underlying assumption 
is that the actual Hamiltonian is complicated in any basis excluding some exceptional ones 
which form a manifold of measure zero. This led, for a time-reversal invariant physical sys- 
tem, to the Gaussian Orthogonal Ensemble (GOE) and, for different global symmetries, to 
other "canonical" ensembles. 

In a sense, canonical Gaussian ensembles describe the extreme degree of universal chaotic- 
ity, and the real interest is in understanding when, and to what extent, this limit is realized 
in actual physical systems. The spectra of quantum systems whose classical counterparts 
are chaotic are well modeled by the GOE: they have the same spectral fluctuation properties 
. This leads us to accept RMT as a working definition of quantum chaos: a quan- 
tum system is deemed chaotic if its spectra exhibit the same local fluctuation properties as 
those of the appropriate Gaussian ensemble. This definition frees us from the need to work 
backwards from classical mechanics. 

The correspondence between the GOE and nuclear spectra has been verified many times, 
over a surprising range of energies. Careful analysis showed that RMT agreed well with 



the neutron 



0, y, Q, 



ll| and proton [1^ resonance data for various isotopes. The range 



of energies at which the nucleus exhibited signatures of chaos was extended all the way to 
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the ground state region in odd-odd nuclei and to two-quasiparticle threshold in even-even 



nuclei [13, 14, 15, 16|. Furthermore, shell model calculations exhibit many of the fluctuation 

nn n 

properties of the GOE |17|, 118|; for an account of tests of RMT in nuclei see [19[. Ongoing 
experimental high-precision studies of neutron resonances in heavy nuclei [2^ require more 
attention to the details and statistical justification of RMT in its practical applications. 
Thus RMT provides us with an arsenal of tools, in the form of certain useful statistics, to 



analyze neutron resonance data 



21|, the largest available body of nuclear spectra. A statistic 



is a number, W, which can be computed from a sequence of levels, and whose mean, (W), 
and variance, Var(W) = {W'^) — (W)'^, are calculable from theory, and which has a small 
deviation from the theoretical value when the theoretical model is a valid one. In our case, 
the theoretical model is that the local fluctuations of the energy levels of excited nuclei are 
described by RMT. To use these statistics, the experimental energies must be rescaled to 
the uniform level density. This process, called unfolding, removes secular variations of the 
spectra, leaving just the fluctuations about the mean values, thus allowing spectra from 
different physical systems to be compared. 

Among various statistical measures of spectral fluctuations, the neighboring level spacing 
distribution P{s), and the Dyson-Mehta A^lL) statistic that quantifies the fluctuations of 
the number of levels in a given spectral interval, are the most useful in practice. Even their 
qualitative features allow one to quickly get a first glimpse of the character of underlying 
dynamics, — regular, chaotic or intermediate. The quantitative analysis provides more 
detailed characteristics. In order to apply such measures successfully, one needs to know 
the completeness of a fragment of an empirical spectrum and its purity. Missing levels or 
contamination by levels of different symmetry classes leads to distortions of the statistics. 
Apart from that, an important question is that of the ergodicity of those statistical measures. 
The predictions of RMT refer, as a rule, to the ensemble average of the quantity of interest. 
Experiment, on the other hand, typically gives the spectral average of the quantity inside 
an observed subset of the large, formally infinite, spectrum, which is taken as a random 
representative of an ensemble. The ergodic property means that the same results are valid 
for different fragments of a given spectrum as well as for the average over the ensemble. 

The Dyson-Mehta A3(L) statistic will be the focus of this work. We study fluctuations 
of this statistic, its ergodic properties and sensitivity to missing levels and impurities. The 
point of contact with experimental data will be neutron resonance data. There are many 
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other systems that lend themselves to an RMT analysis. The sizes of typical data sets varj 
from system to system. The spectrum of electromechanical vibrations in a quartz block 
for example, can have many hundreds of levels, as can data from superconducting microwave 
cavities. Many of the calculations in this paper have in the hundreds. The main results 
are still applicable to neutron resonance data, even though the majority of experimental data 
have less than 100 levels. We provide an error estimate on the calculation of the fraction of 
missed levels for data sets with = 100. 

In Sec. iniwe define the A3(L) and discuss its ensemble average. In Sec. [Tlllwe discuss 
spectral and ensemble averaging, the ergodicity of the A3(L) statistic, and the correspond- 
ing uncertainties. In Sec. IIVI the calculation of GOE spectra, A3(L), and the unfolding 
procedure will be described. A3(L) is calculated exactly, with no numerical minimization 
procedures. Following that Sec. |V] deals with the calculation of the uncertainties. An anal- 
ysis of actual neutron resonance data with the maximum likelihood method is described in 
Sec. I VII In Sect. IVIII we perform a A3(L) analysis of neutron resonance data and compare 
it with the maximum likelihood method. We summarize all our findings in the Conclusion. 



II. A-i{L) STATISTIC: DEFINITION AND ENSEMBLE AVERAGE 

By definition, 

A3(L) = {mmA,B { /f dE' [ M{E') - AE' - Bf ) (1) 

= (A^^)) , (2) 

where we use the notation (x) for the spectral average of x. This is a measure of the average 
deviation of the spectrum on a given length L from a regular "picket fence" spectrum of a 
harmonic oscillator. M{E) is the cumulative level number (the number of levels with energy 
less than or equal to E). The angle brackets in Eq. ([2]) imply averaging over all values of 
i, the location of the window of L levels within the spectrum. A and B are chosen so as to 
minimize A3(L); they are recalculated for each value of i, the starting point of the fragment 
sliding along the spectrum. 

A series of evenly spaced levels would make M{E) a regular staircase, then A3(L) = 
1/12. At the other extreme, a classically regular system will lead to a quantum mechanical 
spectrum with no level repulsion, the fluctuations will be far greater, and A3(L) = L/15. 
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One can also introduce another useful statistic, the level number variance, [2^. It is 

the variance in the number of levels found in an interval of length L. After unfolding the 



spectrum, one expects there to be L± ■^Y?'{L) levels in the interval. For a regular spectrum 
one has T?{L) = L, while for a harmonic oscillator spectrum it is zero. The relationship 



between S^(L) and A3(L) is given in 221] as 

AsiL) = /^(L^ - 2LV + r^)J:^{r)dr. (3) 
Jo 

The asymptotic RMT result for the Gaussian Orthogonal Ensemble (GOE) is 

(4) 



ML) = ^ 



^2 



5 vr^' 

log(27rL)+7--- — 



with 7 being Euler's constant. We stress that this is the RMT value for the ensemble average 
of Al{L), not a spectral average. Putting in values we get A3(L) = (logL — 0.0678)/7r^. For 
the GOE this statistic increases very slowly with L, the levels are crystalized into a rigid 
structure, hence the alternative name "spectral rigidity" for this statistic. 

In this paper we are concerned with detailed properties of the A3(L) statistic and its use 
for the RMT analysis of neutron resonance data. The main issues we will address are the 
fraction, x, of missing levels, and the uncertainties on the A3(L) calculation. A3(L) has been 



applied in this context before. In 23| Monte Carlo calculations were used to see the effect 
of missing levels on pure and mixed GOE spectra. They give empirical graphs that can be 
used to get x, given A3(L) of a specific experimental spectrum. They also give an empirical 
expression for the uncertainties in A3(L) that include effects of both sample size and x. In 
5l an expression for the uncertainties is suggested, and we verify it here numerically. In 
2J] the effect of sample size on A3(L) was examined, and the level spacing distribution was 
deemed a more useful statistic. We reexamine the question here of how to compare A3(L) 
calculated from a set of neutron resonance data with RMT. We will give an exact method 
of calculating A3(L) for an unfolded spectrum, and a consistent approach to comparing the 
experimental result with the theoretical model. 

To illustrate the problem, see Fig. [H where A3(L) is shown for 20 out of a set of 50 GOE 
spectra of random GOE matrices of dimension = 4000. Notice the lines tend to be quite 
smooth, but there is a considerable spread. The mean value of A3(L) for these 50 spectra 
is shown in green in Fig. [H but is not visible as it lies on the theoretical curve, Eq. (jlj). 



This average of the red lines is A3(L), the ensemble average of A3(L). In this way we have 
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FIG. 1: (Color online) Calculated A3(L)-statistic for an ensemble of 50 GOE spectra with = 
4000, results for 20 spectra are shown as thin lines (red). The ensemble average, A3(L), for all the 
50 spectra lies on the theoretical thick (black) curve, differing from it by ~ 2 x 10"'^ over the whole 
range of L. 

recovered the theoretical result, the difference between theory and calculation ^ 2 x 10^'^ for 
the large range of L shown. The source of the discrepancy for a specific spectrum is simply 
the natural spread in A3(L) values for different spectra. This begs the question: if a A3(L) 
calculation on some experimental data gave one of the lines in Fig. [T]what conclusions could 
be drawn about the purity of the spectra, missing levels, etc. We need an "uncertainty" 
to define a confidence interval centered on the theoretical line, so we can make meaningful 
comparisons with the data. 
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III. ERGODICITY OF As{L) 



The validity of a comparison between the spectral average of a quantity with the theoret- 
ical ensemble average depends on the quantity being ergodic. For a clear and more detailed 



a. 



account of this topic see 

Consider an observable X{E), which is some function of energy. In RMT, this ob- 
servable would be calculated by evaluating X{E) for fixed E and averaging over an en- 



semble of spectra to get the ensemble average, X{E). The variance of X{E) is written 
Vare(X) = X{Ey — X{E) ^. We will write the standard deviation as aeX, or simply cJe, 



the subscript indicating ensemble averaging. Dyson 211] derived the variance of A3 in GOE 
to be Vare(A) = 1.169/7r^ = 0.110^, so de = 0.110. On the other hand, experimentally 
one is dealing with an interval of the spectrum over an energy range (determined by the 
experiment) [E,E + AE], and one calculates the spectral average within that range, as a 
running average over the energy, 

1 rE+AE 

{X{E)) = ^J^ X{E')dE'. (5) 



Ergodicity is equivalent to the statement {X{E)) = X{E). 

Take, for example, the nearest neighbor level spacing s = -Ei+i — Ei. It is easy to verify 
that s is ergodic, see Fig. [21 In this case the probability density for s is the same within a 
spectrum as it is in the ensemble. The formal requirement for X{E) to be ergodic is that 



VaTe{X{E)) = {X{E))^ - {XiE))"^ as AE ^ 00. As it stands it is not particularly 
useful, what we need is the behavior of VaTe{X{E)) for finite data sets. Specifically, we need 
an expression for the uncertainty in the quantity after replacing ensemble averaging with 
spectral averaging. We will refer to this quantity as a, with Varf.{X{E)) = a^. 

In application to A3(L), the energy dependence is in which indicates the location Ei 
of the window of L levels. The spread in the individual lines in Fig. [1] is a, while their 



average is (A3(L)), or A3(L), which is the notation we will use. It is cr that will determine 
the sensitivity of A3 as a tool for detecting missing levels. In a calculation of A3(L) on a 
spectrum of levels, we take a spectral average of A\{L), where the average is taken over 
the N — L possible locations E^ of the window of L levels. Brody et al. [5] discuss the 
situation where a quantity X is calculated over p non-overlapping intervals. They suggest 
that Vare(X)p = VaTe{X)/p. They call this the Poisson estimate. In the case of the A3 
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FIG. 2: (Color online) The distribution of level spacings, s, taken from one = 6000 GOE 
spectrum, left, and the distribution of the center level spacing, E^qq — £299, taken from an ensemble 
of 6000 GOE spectra with N = 600. The distributions are the same. The statistic s is ergodic. 

statistic, there are N/L non-overlapping intervals available for each L, and the Poisson 



estimate would prescribe an uncertainty of a = O.lly L/N. We will verify this numerically 
for the GOE. In practice the intervals used in the calculation overlap, as we allow i to take 
all available N — L values. The values of A3(L) are highly correlated in this case however, 
and the Poisson estimate is still good. It is clear from this result that A3 is a more sensitive 
statistic for small values of L/N. 

In the case of m independent spectra, superimposed in proportions /i, . . . fm, and 



letting A^rn{L) be the spectral rigidity of the m^^ sub-spectrum, we have 21 1 A3(L) 
YliLi ^3m{fiL) . If the m spectra are all from the GOE, the ensemble variance is cXg 



(O.llOm)^, and the Poisson estimate then gives a = 0.110m y L/N. Based on this number 
A3(L) can be used to distinguish between spectra with m = 1 and m = 2 independent 
sequences present. However this statistic is not sensitive to the actual mixing fractions. We 
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FIG. 3: (Color online) The dotted (blue) curve is A^L), the ensemble average of AKL), for 



indistinguishable from the RMT prediction, solid (black) line. The lower (red) curve is o"e, the 
standard deviation of A^L). It agrees with the RMT result of 0.11. 

have assumed here that the proportions /i, /2, ■ ■ ■ fm are independent of energy. This may 
not be the case, in neutron resonances, the fraction of intruder p-wave resonances may grow 
with energy. 

To calculate the ensemble average, A^L), 1000 matrices of dimension 3000 were made. 
Each of these spectra was unfolded. The value of z = y — was chosen to locate the window 
of L levels squarely in the middle of the spectrum, A3(L) was calculated, and the resulting 
A3(L) matched the RMT result, see Fig. [31 As a double check, the process was repeated 
for the value of i = ~ f ' ^hich locates the window of L levels on the left edge of the 
spectrum. A discrepancy here would undermine our unfolding procedure, see Sec. HVj but 
the results were identical. 

IV. CALCULATION OF A^iL) 

To realize the GOE, we generated random matrices with normally distributed matrix 
elements, and 



i = 1200 - ^; the ensemble of 1000 spectra with N = 2400 was used. The ensemble average is 





9 



for the off-diagonal and diagonal elements, respectively. We chose cr = 1. Each of the 
matrices has an approximately semicircular level density, with p{E) = \/AN — E'^, for l-E] < 
2\/iV, otherwise (actually there are deviations from the semicircle at the edges, see Mehta 



25|) 



In order to apply the results of RMT to real data, the empirical spectrum must be 
unfolded separating the fluctuations of the spectra from secular behavior and expressing 
all energies in units of the local average level spacing. This is achieved by first extracting 
the cumulative level density ^/{E), which will be a staircase function, from the raw data, 
and fitting it to a smooth function, C,{E), either numerically, or analytically. Then, using 
this function, the j^^ level of the unfolded spectrum is simply ^{Ej). The resulting unfolded 
spectrum has a uniform level density p{E) = 1 so that data from high density regions can 
be compared with data from low density regions. Furthermore, spectra from very different 
physical systems can be compared. 

Integrating the semicircular GOE level density gives 

m =j^^p{E')dE' 

= -.^ {vm^w) + h eVan^ + IN. (7) 

To decide between using the analytical form, or a numerical fit, the curve fitting tool in 
MATLAB was used to get the best values for a, b, and c in the parametrization 

^(E) = a-N tan-^ ( ]+b —EVAN - E^ + c -N. (8) 

The result is, with 95% confidence bounds, a = 0.9992 (0.9992,0.9993), &=!(!, 1.001),and 
c = 1 (1,1). Using these values for a test spectrum with = 500 made a difference in 
A3(L) of 10^^ at L = 250, so the theoretical (semicircle) result, Eq. ([7]), was used in all 
the calculations of A3(L) for GOE spectra, without any fitting. To realize ensembles of 
GOE spectra of various spectrum size, A^, we first made an ensemble of 3000 unfolded GOE 
spectra with dimension 3000. Getting a spectrum of size A^ is now a matter of taking the 
middle A^ eigenvalues from one of these. In what follows, all calculations are performed on 
unfolded spectra. 

We rewrite for convenience the definition of A3(L) for a spectrum: 

rEi+L 

miiA R — / 



A3(L) = (mmA,B ^ J^^^^ dE' [ M{E') - AE' - Bf \ . (9) 
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The integration is over a window of the spectrum, of length L levels, starting at energy Ei, 
the i*^ level. Note that some authors choose the limits of integration to be E-l, E-,l. The 

2 2 

integral is evaluated for every starting energy E^, with 1 < i < N — L. A and B are chosen 
so as to minimize the integral for each position of the window of L levels, i.e. for each value 
of i. The precise values of A and B will be given in terms of Ei, no numerical minimization 
is necessary. 

Substituting Mi^E) = i, E^ < E < into ([9]), our job is reduced to finding the mean 

of the quantity Al{L). In evaluating 

AUL) = I /^'"^^ dE' [ M{E') -AE'-B]\ (10) 

L JEi 

consider the integral between two adjacent levels, 
Ai(L) = y r^^' dE'{j-AE'-By = ]rx{C+VA^+WA+XAB+YB+ZB^), (11) 

L/ ■ ■ J Ei Li 

where 



C = j\E,^,-E,„ 
j=i 

1 

j=i 

X = (-^»>L - Ei), 

Y = ^y'-2j{E,+^-E,), 
j=i 

Z = (Ei^L — Ei). 

Using the constraints d{A\)/dA = and d{A\)/dB = 0, we come to the following 
expressions for A and B that minimize A\{L): 

XY-2WZ WX-2VY 

AVZ-X^ ' AVZ-X^ ■ ^ ^ 

Given an unfolded spectrum, A3(L) can be calculated exactly. Bohigas and Giannoni derived 



V = ^{eIl-e! 



a similar expression for the case of L = — 1 in [26j. It is interesting to note that many 
early investigations of A3(L) deal with just this special case, and hence distinction between 
a spectral average and an ensemble average didn't arise, there being only one value of the 
statistic per spectrum. 
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V. CALCULATION OF a{L) 



To calculate a, the uncertainty for A3(L) of a specific GOE spectrum, we take 1000 GOE 
spectra, each with levels, and calculate A3(L) for each, a is the standard deviation of 
those 1000 numbers for each L, and this is what we should use as the uncertainty for A3(L) 
for that value of L. In Fig. IHwe show the A3(L) and cr for a range of A^. The most important 
feature is that a < for small values of L. When L = N then we have one number per 
spectrum, and a = cTg as expected. 

To compare different size spectra, we plot a vs. L/N, in Fig. [51 A numerical fit to the 
range < L/N < 0.2 gives 

a = 0.1126(L/A^)°-^°°3 - 0.0045 ^ Oe^L/N. (13) 

This is a verification of the Poisson estimate. The dashed line in Fig. [5]is a = af.\J L/N . 
The Poisson estimate was stated for the case of non-overlapping intervals. In Eq. (10) this 
would mean i = l,L + l, 2L + 1,.... In our calculations all values of i were used, this 
is equivalent to doing L sets of calculations for non-overlapping intervals, with the initial 
values of i in each set going from 1 to L, and then taking the average of these L numbers, 
so the Poisson estimate is still valid. Note that these "uncertainties" are not to be confused 
with any experimental uncertainty or computational issue. Given an unfolded data set, the 
A3(L) statistic is calculated exactly, as was shown in Sec. HVl 

For the case of a mixture of m = 2 independent GOE spectra, the ensemble result is 
(Te = 0.22. A calculation of a for 500 spectra each with 500 levels, made by mixing two 
GOE spectra, with /i = 0.4, gives a = 0.239(L/A^)°-^22 _ o.oi4, which is consistent with 
a = 2 ae\J L/N . In Fig. [6] we show A3(L) with the empirical uncertainties. The statistic 
can clearly distinguish between m = 1 and m = 2 spectra, and it is most sensitive for 
small values of L. A-^^L) is not sensitive to the actual value of /i, and cannot be used 
to distinguish between f\ = 0.25 and /i = 0.4375. These are the values for the fraction 
densities one would expect for neutron resonances on target nuclei with spin jtar = 1/2 and 
7/2, based on the 2j + 1 degeneracies of the resulting compound nucleus spins jtar ±1/2. The 
statistic is sensitive to p-wave neutrons in the beam being mislabeled as s-wave for spin-1/2 
targets, as this would introduce multiple sequences of levels into the data set. 

To determine the completeness of a spectrum with A3(L), we need to know its behavior 
for GOE spectra with a fraction of levels, x, depleted. We denote this by A3(L;a;), and 
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FIG. 4: (Color online) The ensemble average, A3(L) vs. L is plotted for various values of N. Each 
ensemble had 1000 matrices. The lower lines are the standard deviations of the corresponding 1000 
values of A3(L). Notice that for the maximum value of L, a = 0.11, which is the ensemble result. 
This is expected, as there is only one such value per spectrum. 
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FIG. 5: (Color online) Here we plot a from Fig. |3]vs. L/N, for various ensembles. The top (black) 
line is cie, in agreement with the RMT result of a constant 0.11 . The dashed line is ae\/L/N. 

it was found empirically for ensembles of 500 spectra with dimension and depletion x. 
The dimension, A^, refers here to the number of levels after the fraction x was randomly 
removed. The ensemble average of A^{L;x) and its standard deviation a{N,L;x) were 
calculated. The results are shown in Fig. [7] for N = 200, with uncertainties. We see 
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FIG. 6: (Color online) A^lL) is plotted with uncertainties for the cases of an ensemble of 500 pure 
GOE spectra, lower lines (black), and a 500 mixed GOE spectra, upper lines (blue). In the mixed 
case, each spectrum was a superposition of 2 pure GOE spectra in proportions fi = 0.4 , /2 = 0.6. 
In both ensembles, the spectra had N = 500. The uncertainties are empirical, being the standard 
deviation of the 500 values of A3(L) at each value of L. 

from Fig. [S] that the uncertainty here has the form a^iN, L; x) = f{L)/\/l^. The l/\fN 
dependence of crx{N, L; x) means that the statistic is a less sensitive measure of x for lower 
A^. A fit to f{L) was made for a range of x, to be used in practical analysis of the neutron 
resonance data. 

Now we have the tools necessary to compare a real spectrum with depleted GOE spectra 
in a meaningful way. The best value for x from the data will be the one that minimizes 



For practical purposes we need an estimate of the error in x from this method. To this 
end, we calculated the average value of x, and its standard deviation, for 1500 GOE spectra 
with = 100 and 300 levels. The sets were made by randomly deleting 3 levels from a 
spectra of 103 levels, 5 out of 105, 9 out of 109, and 11 out of 111, to get spectra with 
X = 2.91%, 4.76%, 8.26%, and 9.9% respectively. The results are in Table [H Our method 
gave good agreement for the value of x, but the uncertainties were of the same order as x, 
for example, with x = 4.76%, we get a value of (4.89 ± 2.81)%. However, when we tried the 
method for A^ = 300, the uncertainty dropped by a factor of nearly to 1.66%, which is 



[A,{L)-As{L;x)Y 
a{N,L]xf 



(14) 
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TABLE I: The errors in x, the fraction of missed levels, using the A3 statistic. The tests were run 
for 1500 depleted spectra of size N = 100 and 300. The mean value x, and the standard deviation, 
a are given. 

X x{N = 100) cj(A^ = 100) x(iV = 300) a (iV = 300) 

2.91% 3.14% 2.50% 2.71% 1.57% 

4.76% 4.89% 2.81% 5.00% 1.66% 

8.26% 7.95% 3.13% 7.94% 1.74% 

9.91% 9.90% 3.06% 10.0% 1.77% 
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FIG. 7: (Color online) A^(L) vs. L is plotted for depleted spectra size, = 200, with uncertainties. 
Lower line has x = 0% depletion, middle line has x = 5%, and the upper line has x = 10%. 

to be expected. 

An experimental spectrum can be contaminated by intruder levels. In the case of s-wave 
neutron resonance data, this could be from the capture of p-wave neutrons. To examine 
the effect of contamination on A3(L), ensembles of 1000 unfolded GOE spectra with A^ = 
200(1 — x) were made. Each spectrum was stretched by a factor of 1/(1 — x) and then 200a; 
random numbers, uniformly distributed on the interval [0,200], were added. The resulting 
spectrum had a uniform level density, p{E) = 1, a size A^ = 200, and a level of contamination 
of x%. In Fig. Owe see the results for x = 2%, 5%, and 10% (solid lines). The results agree 
well with the RMT prediction of A3(L) = (1 - x)[(logL - 0.0678)/7r2] + xL/15 (dashed 



— 0% depleted 
—5% depleted 
—10% depleted 




15 




200 400 600 800 

L 



FIG. 8: (Color online) \n^a vs. L is plotted for depleted spectra of dimension A^. The values of 
N are 200, 400, 600, and 800, L goes out to its maximum value of — 1. 
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FIG. 9: Here A3(L) is shown for contaminated GOE spectra, upper solid lines. The 

level of contamination, x, is labeled. Upper dashed lines are the RMT predictions. The 

dashed lines at the bottom of the plot are the standard deviations, a, of A3(L) for x = 
2%, the lowest of the three, to 10%, the highest of the three. 

lines). 

We have discussed the spectral average of A3(L), but what about its spread within a 
spectrum? We call the standard deviation of A3(L) for a given spectrum cr^. In Fig. [TO] we 
see the average of ag as a function of L for different spectra sizes, A^. It is immediately clear 
that as is less than ae- Within one spectrum, we expect a smaller spread in the values of 
Al{L) because close values in i mean the windows of L levels overlap, the corresponding 
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FIG. 10: (Color online) The standard deviation of A\{L) within a spectrum, ag, vs. L is plotted 
for various values of A^. The linear parts of the plot are well described by as = CFe{N — L) 

values of A3(L) would be correlated, and (Xs should get smaller as L — A^. A plot of 
(Js versus L/N strongly suggests the falloff is linear, with o"s = cre(A^ — L), see Fig. [TTl 
It is not obvious that the correlations between overlapping windows of levels would give 
this linear behavior. To gain insight into the correlation between A3(L) and A3'^^(L), we 
examined the ensemble average of the square of the difference between these two quantities, 
specifically, 5A{L) = \J[A'^^^(L) — A^i^L)]"^, the variation of A3(L) with respect to i. It is 
expected that this quantity should be closely related to A3(L), and it certainly decreases 
rapidly, SA{L) oc A^{L)/L. The constant of proportionality is = 2.3303, where 

the average was taken over all values of L, see Fig. [121 

It is interesting to see the behavior of as{N, L; x), the spectral average of the variation of 
A3(L; x). It has a strong A^- and x-dependence, and the simple rule y/N — L of the x = 
case is lost. In Fig. [13] we show that the A^-dependence of crs{N, L; x) for 5% depletion peaks 
at L ^ A^/2. Strong correlation in overlapping windows L levels wide would reduce this 
number as L increases, and in Fig. [TH the results for A^ = 200 and x = 0%, 2%, 5%, and 
10% show that the maximum moves to higher values of L as x increases. 
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FIG. 11: (Color online) as vs. L/N is plotted for various values of N. The linear parts of the plot 
are well described by 0"<j = cJe(l — -L/A^) 




FIG. 12: (Color online) The quantity 5A{L), as defined in the text, is plotted on a log-linear graph. 
It is well approximated by 2.3303A3(L)/L (smooth line). 

VI. LEVEL SPACING ANALYSIS 

The nearest level spacing distribution can be used to test for missing levels, and mixing 

n 

of different sequences. Here we will follow the work of Agvaanluvsan et al. [27|] , where the 
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FIG. 13: (Color online) The rescaled standard deviations for systems of various sizes, as{N, L;x), 
with 5% of levels depleted. Each line represents the average for 500 systems of a particular size. 




L/200 

FIG. 14: (Color online) The rescaled standard deviations for N = 200, as{N, L; x), with varying 
degrees of depletion. Each line represents the average of 500 systems. 

maximum likelihood method is used to find the fraction of missing levels in a sequence. We 
will test the method on GOE spectra depleted by hand. Single sequence neutron resonance 
data will then be analyzed, and the results compared with those of a A3 analysis. 

The level spacing distribution for a complete GOE spectrum is given by the Wigner 
surmise. 

Pis) = \se--^'l\ (15) 
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FIG. 15: (Color online) The maximum likelihood method (MLM) is tested on a superposition of 200 
GOE spectra, each of length = 250. A fraction x of each spectrum was randomly removed, and 
the MLM was used to recover the number. Here are the results. The MLM slightly overestimates 
the fraction missing. The natural log of the likelihood function, ln(£), is plotted against x. 

where s = S/ D, S being the spacing between adjacent levels, and D is the average spacing. 
We have unfolded all spectra involved, so D = 1. If a level is missing, then two nearest 
neighbor spacings are unobserved, while one next-to-nearest spacing is included as a nearest 
level spacing, when it should not be. Furthermore, if 1 — a; is the fraction of the spectrum 
that is observed, then Z^obs? the experimental value for the average spacing, is related to the 
true value hj D = {1 — x)Dohs- Agvaanluvsan et al. show that 

oo 

P(s) = ^(1 -x)a;'=P(/c;s), (16) 

fc=0 

where P{k; s) is the distribution function for the k^^ nearest neighbor spacing, Ek+i — Ei] 
for = this reduces to the Wigner distribution, P(0; s) = P{s). 

The maximum likelihood method will be used to find the best value for x that maximizes 
the likelihood function C = Y\iP{si); the product is over all the observed spacings. In 
practice it is easier to maximize ln(£) = X]iln-P(sj). The functions P{k; s) are complicated 
to derive. For k ^ 0, the functions were fitted to the empirical distributions from the 
superposition of 5000 GOE spectra, each of length = 5000. The procedure was tested on 
200 GOE spectra, each of length N = 250, with 3%, 5%, 10% and 13% of levels randomly 
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TABLE II: The errors in x, the fraction of missed levels, using the maximum likelihood method. 
The tests were run for 1500 depleted spectra of size = 100 and 300. The mean value x, and the 
standard deviation, a are given. 

X x{N = 100) cj(A^ = 100) x(iV = 300) a (A = 300) 

2.91% 3.32% 2.42% 3.44% 1.49% 

4.76% 5.02% 2.56% 5.29% 1.43% 

8.26% 8.45% 2.48% 8.74% 1.32% 

9.91% 10.09% 2.33% 10.32% 1. 




0.05 0.1 0.15 0.2 

X 



FIG. 16: (Color online) MLM results are shown for some representative data sets. The indication 
is that the ^'^^U data is incomplete, with 5% of levels missing, but the subset of the data consisting 
of the first 69 levels looks complete. There appears to be 3% of levels missing from the ^'^^W data, 
with 0% missed from the ^'^Ni data. 

removed. The results are shown in Fig. [151 The method systematically overestimated x by 
about 0.5%. 

In order to get an estimate of the error in x from this method, the procedure of the 
previous section was repeated, and the results are shown in Table [Tll The agreement with 
the A3(L) method is encouraging. The A^lL) method seems to be slightly more accurate, 
but the uncertainty in x is slightly smaller for the MLM. 

Next the procedure was applied to the real neutron resonance data. The data sets are 
described in the next section. In Fig. [16] we see some typical results of the MLM analysis. 
For the ^'^^U data, the full set of 81 levels looked incomplete, to the tune of about 5%, while 
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the first 69 levels looked complete. The ^^^W data looked like there were 3% of the 68 levels 
missing, while the 63 levels in the ^^Ni data looked like a complete set. 



VII. NEUTRON RESONANCE DATA 



The neutron resonance data for a wide range of isotopes are now widely available |29 |. 
Data sets with spin-0 target nuclei were chosen for analysis. They afforded the simplicity of 
not having the fractional densities /i and /2 of a superposition of two independent spectra, 
as all levels have spin 1/2. There is a conventional assumption here that the capture is 
dominated by s-wave neutrons, so no higher angular momentum resonances intrude. Eleven 
isotopes were analyzed, and the data sets are described in Table IIIII 

The cumulative level number gives the first indication of the purity of the data. Kinks in 
N'{E) leading to smaller slopes would suggest a section of data where levels were missing. 
Using this as a guideline, some data sets were split into subsets. For example, in Fig. [T71 
there is a kink in J\f{E) at the 70*^ level, so we analyzed the subset of levels Ei E-^q. The 
percentage of missing levels was estimated with the maximum likelihood method. A3(L) 
was then calculated for each set, and compared with the MLM results. In Fig. [17] the three 
quantities are shown for two representative isotopes, ^^^Gd (blue), and ^^^Gd (red). The 
results are summarized in Table UTTl The agreement between the two methods is promising. 

There is a range of behaviors exhibited by A3(L). It can be flatter than the RMT value, 
which suggests an artificially rigid spectrum. This type of behavior is typical when the 
unfolding procedure employs too specific a function for ^/{E), and all the fluctuations are 
washed out, leaving the unfolded spectrum too rigid. It can grow too rapidly, which could 
mean that there are more independent spectra present than thought. This would happen if, 
for example, p-wave neutrons were contributing to the data. 

In the presence of sequences of energy levels with different spin labels, we would have a 
superposition of independent spectra. This is the case, for example, when s-wave neutrons 
are incident on a target nucleus with spin jtar, so that the resulting resonances have spin j± = 
jtar =t |. The level density is relatively constant over the small energy range of the data, and 
one would naively expect that the subspectra have fractional densities /i = _^^j~j^2j++i) 
and /2 = 1— /i- For m independent spectra, with fractional densities /i , /2 • • • /m, and letting 
A^rniL) be the spectral rigidity of the m*^ sub-spectrum, we have A3(L) = X^Ili ^smifiL). 
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FIG. 17: (Color online) A3(L) vs. L for the experimental data. The black lines are the GOE 
values for 0%, 2%, 5%, and 10%, starting from the lowest curve. The ^^^Gd data is in blue, and 
the ^^^Gd data is in red. 



A preliminary analysis of systems of mixed independent spectra was done on the ^ss-q 



data [28(1 . The s- resonances have spins j = 3 and 4, making the data a superposition 
of two independent spectra, with fi = 0.4325. This data set exhibited all the behaviors 
described above. This data set, the biggest by far, was split into four sections corresponded 
to linear regions on the M{E) curve. The data was unfolded by fitting sections of N'{E) to 
a straight line. The first section had 950 levels and the energy range from to 510 eV. It 
exhibited excellent agreement with the GOE result for mixed (/i = 0.4) spectra with 4% 
depletion. The next sets were levels 1050 to 1450, with range 578 eV 920 eV, 1700 to 
2700 with range 1118 eV 1995 eV, and the fourth set had levels 2750 to 3150, with range 
2022 eV 2240 eV. The results are summarized in Fig. [191 Although we do not have a 
MLM comparison for the case of two independent spectra, the /S.^{L) analysis suggests that 
4% of the levels were missing in the first 960 levels. 



VIII. CONCLUSION 



We have reexamined the possibility of using the Dyson-Mehta A3(L) statistic as a di- 
agnostic for the spectra of quantum systems in the chaotic regime, where the fluctuations 
can be modeled by Random Matrix Theory. Originally it was regarded as a tool of much 
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TABLE III: The results for x, the percent of missing levels in the data. 



Isotope 


MLM 


ML) 


N (# levels) 


subset 




2% 


4% 


64 


1 


36 


50Cr 


0% 


8% 


64 


37 - 


64 




3% 


0% 


63 


1 


50 


58N 


0% 


0% 


63 


All 




i52Sm 


3% 


2% 


91 


1 


70 




8% 


12% 


128 


1 


70 




4% 


2% 


161 


All 




158Gd 


11% 


7% 


93 


All 




158Gd 


0% 


1% 


93 


1 


60 


158Gd 


12% 


15% 


93 


61 - 


93 


182^ 


3% 


2% 


68 


All 




234u 


6% 


4% 


118 


1 


75 


234u 


7% 


3% 


118 


76 - 


118 


234u 


9% 


8% 


118 


All 




236 u 


0% 


1% 


81 


1 


69 


236 u 


5% 


8% 


81 


All 




240p^ 


0% 


3% 


267 


1 


100 


240pu 


8% 


12% 


267 


171 


^ 267 


242pu 


11% 


13% 


67 


All 





more limited resolution, due to the large variance [2l|. This is not the full picture, and 
the definition of the statistic and the original results are reinterpreted. An examination of 
the ergodicity of the statistic lead to the spectral properties being separated from ensemble 
properties. We see that it is inappropriate to use the large and constant variance of the 
statistic as an uncertainty. Uncertainties in A3(L) were found empirically for pure and de- 
pleted GOE spectra. In doing this, the Poisson estimate of Brody et al. was verified. The 
A^iL) statistic was used to determine the percentage of missing levels in neutron resonance 
data. The results were compared with the maximum likelihood method of Agvaanluvsan et 
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FIG. 18: (Color online) A^iL) vs. L is plotted for the first 960 levels of the ^^^V data. The 
uncertainties are empirical for mixed (/i = 0.44) GOE spectra, with 0%, 4%, and 10$ depleted. 
The data is consistent with there being 4% of levels missing. 




FIG. 19: A3(L) vs. L is plotted for the four subsets of the ^^^U data.The plot suggests that 
the middle sections of the data are less pure than the first set, but the highest energy set looks 
artificially rigid, like a harmonic oscillator spectrum. 

al and the agreement was good. The method was apphed to ^^^U data with curious 
results. Various sections of the data displayed different behavior, raising questions about 
the purity of the data, the accuracy of the angular momentum assignment, and the appro- 
priateness of the modeling the data with RMT. The behavior of ag, the spectral spread of 
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A\{L), was examined, and new properties were described which need a more detailed study. 
It is hoped that the A3(L) statistic will be considered as a useful diagnostic in the RMT 
arsenal for data analysis. 
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